The 11^'^ Asian-Pacific Regional lAU Meeting 2011 

NARIT Conference Series, Vol. 1, ©2012 

S. Komonjinda, Y. Kovalev, and D. Ruffolo, eds. 



Of Winds and Waves 



John G. Kirk^ and loanna Arka^ 

^Max-Planck-Institut fiir Kernphysik, Postfach 10 39 80, 69029 Heidelberg, 
Germany 

^Institute de Planetologie et d'Astrophysique, University of Grenoble, 
Grenoble, France 

E-mail: john.kirk@mpi-hd.mpg.de 

Abstract. Recent work on the properties of super luminal waves in pulsar 
winds is summarized. It is speculated that these waves play an important role 
in the termination shock that divides the wind from the surrounding nebula. 



1. Introduction 

In addition to pulses of radiation, rotation-powered pulsars are thought to 
emit a wind that powers the diffuse radiation or pulsar wind nebults (PWN) 
observed around many of them in the radio, optical. X-ray and gamma-ray bands. 
Synchrotron emission and inverse Compton scattering are the mechanisms most 
likely to be at work in PWN, and it has been evident for over half a century that 
not only the energy, but also both the relativistic electrons (or positrons) and 
the magnetic fields in these nebulae must be supplied by the central star[l]. 

The basic idea is that the rotation of the neutron star couples to the wind 
through the strong magnetic field anchored in its crust [21 [3]. As well as 
electromagnetic fields oscillating at the rotation frequency of the pulsar, the wind 
contains a DC (phase-averaged) component of the magnetic field, electrons and 
positrons created by cascades in the magnetosphere and, perhaps, a relatively 
small number of electrons and ions extracted from the stellar surface. As it 
propagates away from the star, the ram pressure of the wind decreases, and, 
roughly where it equals the ambient pressure, a termination shock is formed. 
Rees & Gunn [3] simply assumed that the waves are absorbed at the termination 
shock, leaving behind the energized particles and magnetic fields that fill the 
nebula. Kundt & Krotscheck [3], however, noted that electromagnetic waves 
could also be reflected by the shock and build up inside it. This raises an 
important question concerning the structure of the termination shock and the 
wind it encloses. Because the waves have large amplitude, they interact strongly 
with each other and with the particle component. Rather than just a linear 
superposition of DC fields, particles, outward propagating waves and reflected 
waves, one must, therefore, look for self-consistent solutions, containing all of 
these, that match the outer boundary conditions. 



2. The a problem 

Estimates of the energy density in magnetic field at the light cylinder tl = c/uj 
(with uj the angular frequency of the pulsar) suggest that the wind is energetically 
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dominated at launch by the waves and/or DC fields. On the other hand, outside 
the termination shock the particle pressure is comparable or larger than the 
pressure exerted by the fields (for reviews, see [SI [6]). The implied conversion 
of electromagnetic into kinetic energy is quite natural in many situations, 
for example, when an MHD outflow is collimated and accelerated into a jet. 
However, this kind of conversion does not appear possible in a pulsar wind, 
where the poloidal flux threads the neutron star [6J. In such the MHD 

equations suggest the flow remains radial and its (supermagnetosonic) speed 
stays constant until the termination shock is reached. 

This, then, is the "cr problem" (the ratio of Poynting to kinetic energy flux 
is conventionally denoted by cr): if the pulsar wind can be described by the 
equations of ideal MHD, there is no way in which a large value of a near the 
star can be reduced to unity or below before the termination shock is reached. 
Furthermore, if the termination shock is simply a discontinuity that obeys the 
usual jump conditions in an otherwise ideal MHD flow, the downstream plasma 
remains magnetically dominated with large a. 

3. Dissipation in current sheets 

To solve this problem, we obviously need to go beyond the ideal MHD description 
[7]. One possibility, suggested by Coroniti [8] and Michel [9], is to allow for 
dissipation in a current sheet embedded in a wind that is otherwise a cold, 
MHD flow. This is the "striped wind" picture. Conversion of the energy flux 
from fields to particles by dissipation in the sheets causes the flow to accelerate 
[lUl \TT\ [T^ — in the simplest model the bulk Lorentz factor rises as T oc r^/^. 
Consequently, time-dilation reduces the effective dissipation rate seen in the lab. 
frame. It is possible to place upper and lower limits on the dissipation rate in 
order to constrain the characteristic radius at which the oscillating components 
of the fields are annihilated. The results depend on the pair-loading of the flow, 
described by the parameter 

where L is the luminosity and M the mass-loss rate. For the Crab, dissipation 
occurs before the termination shock is reached only if < 10'^, a value somewhat 
lower than conventional estimates |13j . 

Embedded current sheets also change our picture of the termination shock. 
On the basis of an analytical model and ID PIC simulations, Lyubarsky & 
Liverts [I^ and Petri & Lyubarsky [15] suggested that an MHD shock sweeping 
through the cold portions of the striped wind would drive reconnection in the 
embedded current sheets and dissipate a substantial fraction of the magnetic 
energy, provided the resulting heated electrons were unmagnetized. For cr S> 1, 
this is equivalent to the condition 

> a (2) 



where the strength parameter a of the incoming wave (striped wind) is the ratio 
of the quiver frequency of the electrons eB /mc to the angular frequency to of the 
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wave, (see [T6])I!] Since a decreases as 1/r in spherical geometry, ([2]) amounts to 
a lower limit on the radius at which such a process could operate: 

r > rcrit (3) 
« 3.4 X lO^^rLLgf 

where (L38/47r) x lO^^erg/s is the luminosity per unit solid angle carried by the 
wind. 

On the other hand, Sironi & Spitkovsky [18] recently considered driven 
reconnection in the striped wind using 2D and 3D PIC simulations. Although 
their investigations were confined to the range 0.05 < r /rcrit < 4, they found 
substantial dissipation in all cases. They also noted the appearance of high- 
energy particles accelerated by the first-order Fermi mechanism, and, on the 
basis of a comparison of the nonthermal particle spectrum with the observed 
synchrotron spectrum, suggested that, in the case of the Crab, the termination 
shock should be located where r ~ rcrit/3. If this interpretation is correct, 
the radius of the termination shock inferred from X-ray observations (~ 0.1 pc) 
implies /i ~ 50. This is much smaller than conventionally estimated, implying a 
very high pair loading and a mildly relativistic T ~ 10 outfiow. An approximately 
spherically symmetric wind with these parameters would vastly overpopulate 
the nebula with electrons and positrons. However, if we interpret the wind 
parameters and, therefore, the radius of the termination shock, as latitude 
dependent, it cannot be ruled out that some part of the wind, for example the 
equatorial plane, contains a relatively dense, slow outfiow of the type suggested. 



4. Charge starvation 

As well as giving rise to dissipation in current sheets, non (ideal-)MHD effects 

can also manifest themselves via charge starvation, which can arise, for example, 

when the required currents demand relativistic motion of the current carriers 

[T9l [20] . A relativistic two-fiuid (electron-positron) approach can be used to 

model this effect. In the simplest case, with cold fiuids and charge neutrality, 

both subluminal and superluminal nonlinear waves can be found [16j. 

The subluminal wave resembles the striped wind, except that the current sheet 

^ 2 

is replaced by a static shear, with B constant (i.e., a circularly polarized wave). 
The phase profile of the shear can be chosen arbitrarily, but it is monochromatic 
in the simplest case of constant density. It is interesting to note that, as for the 
striped wind, the non-MHD effects cause also this solution to accelerate with 
radius, in this case with T oc r, despite the fact that no dissipation process is 
involved p!9|. However, at least for the monochromatic wave, acceleration starts 
relatively far from the pulsar, where 



a > a (4) 

^ There are several ways to define this important parameter. In the striped wind, where |-B| is 
phase-independent and a ^ 1, there is no ambiguity. For the superluminal waves discussed in 
section[S]we take a — eE-ma^/mcw, where -Emax is the amplitude of the oscillating electric field. 
A covariant, gauge-independent definition for vacuum waves is given by Heinzl & Ilderton [17] 
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Figure 1. The radial four- velocity of the fluid, the strength parameter a, and 
the refractive index n of circularly polarized superluminal waves close to the 
cut-off frequency for /i = 12, o" = 3. 



corresponding to 

r > Trcrit (5) 

The solutions with superluminal phase speed resemble vacuum electromag- 
netic waves, in the sense that the displacement current plays a crucial role. The 



fields are transverse, and satisfy 



E 



B 



The waves and particles are not 



le subluminal solutions, which makes 



> 

locked together in the way they are in t : 
them more promising candidates for matching to a boundary condition at the 
termination shock. There have been extensive investigations of these waves in 
the literature (see [21] and references therein). They propagate only when the 
density falls below a critical value. In spherical geometry, the corresponding con- 
dition on the radius coincides with ([3]). Because of this, they cannot be launched 
close to the pulsar, and can only be realized if a striped-wind-like (subluminal) 
solution converts into a superluminal one either spontaneously, or because it is 
forced to do so by the boundary condition imposed by the termination shock. 
In this case, the conversion process and subsequent radial evolution and damp- 
ing of the mode should more properly be regarded as an integral part of the 
"termination shock" structure. 
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Figure 2. The refractive indices of linearly polarized superluminal waves close 
to the cut-off frequency for fi = 10100, a = 100 as a function of radius normalized 
to 

'"crit (see Eq. [3]) , for different values of the phase-averaged magnetic field 
5. Superluminal waves 

A first step towards understanding the role of superluminal waves is to determine 
the propagation characteristics of those waves to which a given striped (or other 
subluminal) wind can convert. This amounts to solving the two- fluid analogue 
of the shock jump conditions. Here we present a brief summary of some recent 
work on this topic — a full description can be found in [21j . 

Circularly polarized modes with zero phase- averaged field are the simplest 
to analyze, since the jump conditions can be solved analytically [16]. Figure [1] 
illustrates the properties of these waves for parameters /i = 12, o" = 3 (and, 
therefore, F = 3). Although not in the range expected in pulsar winds, this 
parameter set reveals the mode structure near the cut-off frequency particularly 
clearly. The figure shows that for each value of the frequency uj in the lab. 
frame, measured in units of the proper plasma frequency Up, two solutions exist 
that carry the prescribed particle, energy and momentum fluxes: a strong wave 
through which the plasma streams relatively slowly (dashed line) and a weaker 
one (solid line) through which the plasma streams rapidly (u^ is the four velocity 
along the propagation direction). Transition to the the weak wave involves 
almost complete annihilation of the incoming field energy, leading to Ux ~ /i. 
The refractive indices of these modes coincide and correspond to superluminal 
phase speed n < 1. Because of the finite wave amplitude, the lab. frame frequency 
to! exceeds the proper plasma frequency Wp, which equals the cut-off frequency 
for linear waves. 
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Figure 3. The phase-averaged radial four-speed of the electron and positron 
fluids in the linearly polarized superluminal waves shown in Fig. [2l 

Linearly polarized waves are much more complicated, requiring a numerical 
solution to the jump conditions. Nevertheless, the properties illustrated in Fig. [1] 
apply also to these modes. Two solutions are available for each set of "upstream" 
parameters that characterize the corresponding subluminal solution. This is 
shown in Fig. [21 which differs from Fig. [T] in two ways. Firstly, we show the 
refractive index not as a function of frequency, but of radius in the pulsar wind, 
normalized to rent given in Eq. In this representation, the refractive indices 
of the two solutions no longer coincide, since the proper plasma frequencies 
associated with each solution, differ although the lab. frame densities do not. 
Secondly, we show solutions for four different values of the phase-averaged 
magnetic field, given by the parameter b = (^-B^ / (5^)^'^^. It can be seen 
that stronger DC fields push the cut-off radius further out from the pulsar. 
Furthermore, stronger fields display solutions with negative refractive index, 
corresponding to phase velocities directed radially inwards. Figure [3] shows the 
phase-averaged four-velocity of the plasma stream in the radial direction. Note 
that the weak wave has {ux) ~ /U only for zero averaged field. For finite values 
of the DC field, only a fraction of the incoming Foynting flux is carried by the 
oscillating fields and can be annihilated. 

6. Summary 

Because they are powerful enough to evacuate a relatively large cavity, pulsar 
winds provide an environment in which electromagnetic fields dominate the 
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dynamics, and, especially in the outer regions, non-MHD effects are crucial. 
Within the critical radius given by ([3]), energized electrons and positrons, which 
have Lorentz factors ~ /i, can be regarded as magnetized, since their relativistic 
Larmor radius, fimc^/eB (which in this case equals the effective inertial length 
c/u}p) is smaller than the length scale c/co. Here the superluminal modes cannot 
propagate, because the plasma is sufficiently dense to screen out electromagnetic 
waves reffected from the termination shock. In this region, 2D and 3D PIC 
simulations indicate that dissipation could be driven by a standing shock front 
[I8j . However, in the case of isolated pulsars, rent lies well inside the point of 
pressure balance with the external medium, making it doubtful that such a shock 
could be supported. 

Outside Tcrit energized particles are unmagnetized. Driven reconnection at 
the termination shock may still occur [15j, but superluminal waves may also 
play an important role. In section [5] we have briefly summarized some recent 
work [21] on the properties of these waves for parameters thought appropriate 
in pulsar winds, and sketched their possible role. Although the physics of these 
two scenarios remains to be investigated in detail, it seems likely that there 
will be observable implications. Indeed, current sheets have been proposed as 
the source of both the pulsed high-energy emission [22^ [23] as well as the the 
recently observed flares in very high energy gamma-rays from the Crab (24^ [25] . 
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